Incoherent Pair Tunneling as a Probe of the Cuprate Pseudogap 
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We argue that incoherent pair tunneling in a cuprate superconductor junction with an optimally doped su- 
perconducting and an underdoped normal lead can be used to detect the presence of pairing correlations in the 
pseudogap phase of the underdoped lead. We estimate that the junction characteristics most suitable for studying 
the pair tunneling current are close to recently manufactured cuprate tunneling devices. 



PACS numbers: 74.50.+r, 74.40.+k,74.20.Mn, 74.72.-h 

The pseudogap - a depletion of the single particle spectral 
weight around the Fermi energy - is considered to be one of 
the most convincing manifestations of the unconventional na- 
ture of cuprate superconductivity The pseudogap regime 
sets in as the temperature is lowered below a crossover tem- 
perature T* , and extends over a wide range of temperatures 
in underdoped samples While the pseudogap is clearly 
present in the spin channel [||], optical conductivity data [|j] 
suggest that the same mechanism is responsible for the gap- 
ping of the charge degrees of freedom as well. In addition, 
specific heat data [Q] also provide evidence that a gap opens 
below T* . It has been suggested [^-^|] that precursor super- 
conducting pairing fluctuations may be responsible for these 
phenomena; the observations of a smooth crossover from the 
pseudo to superconducting gap seen in angle-resolved pho- 
toemission ^ and scanning tunneling spectroscopy JlO| ] lend 
support to this idea. There are, however, several other compet- 
ing proposals that do not necessarily involve charge 2e pairing 
p"l|]. It is therefore of interest to find an experiment which 
can provide a direct test of the superconducting precursor sce- 
nario. Here we propose and analyze an experiment involving 



incoherent pair tunneling which provides such a test [ 12 1 



The measurement of the pair susceptibility in the normal 
state of a superconductor, is in principle similar [|l3[| to other - 
say, magnetic - susceptibility measurements: we are interested 
in finding out the linear response of the system to a polarizing 
external field. In the present case the role of the external field 
is played by the rigid pair field of a second superconductor 
below its transition temperature, which couples to the fluctu- 
ating pair field of the normal lead. This coupling leads to an 
observable contribution to the tunneling current - the incoher- 
ent pair tunneling current - provided that the normal state has 
sizeable pairing correlations. 

The basic experimental configuration is illustrated in 
Fig. [I] (a). An I-V measurement is made on a tunnel junction 
formed from an optimally doped cuprate superconductor A 
and a non-optimally doped material B, in a temperature range 
> T > T C B . The c-axis is perpendicular to the A and 



B layers which are separated by an insulating layer. Such a 
structure could be obtained by varying the doping concentra- 
tion of a crystal during a layer-by-layer deposition jl4|]. If the 
B lead is underdoped, as indicated in Fig. [j](b), there will be a 
substantial temperature region above T C B in which B will have 
a pseudogap, while if the normal lead is overdoped (denoted 
by B' in Fig. |j](b)), this pseudogap region will be significantly 
narrower. Now, for > T > Tf we can use the supercon- 
ducting pair field of the optimally doped superconductor to 
directly probe the strength of the pairing fluctuations in B, by 
measuring the incoherent pair tunneling contribution I p (V) 
to the total tunneling current I{V) | |T3| , p3|JTr^ ] . If pseudogap 
behavior is associated with strong precursor superconducting 
pairing, the contribution from the incoherent pair tunneling 
I P {V) should extend over a much wider temperature range 
than for the overdoped B' lead, even if for the two are 
equal. Furthermore, if indeed the pseudogap region is char- 
acterized by precursor pairing, the voltage structure of I P (V) 
provides a measure of the frequency dependence of the imag- 
inary part of the particle-particle t-matrix in the pseudogap 
regime. 

The incoherent pair tunneling contribution to the total tun- 
nel current I(V) is shown diagrammatically in Fig. [j] (c) Jlr^]. 
The quasiparticle states in A are labeled with p and p' and 
those in B with k and k': The p lines correspond to the 
Gor'kov functions F p (iu> n ) of the A lead, while the k lines 
represent the single particle Green's functions of the B lead. 
The dots represent the one-electron tunneling matrix element 
V v ,k and t is the particle-particle t-matrix for material B. It 
should be stressed that, due to the very short c-axis coherence 
length in the cuprates, the pair tunneling takes place between 
the two cuprate layers on either side of the insulating barrier. 
In order to analyze the experimental requirements, we con- 
sider a circular Fermi surface and assume that the pairing in- 
stability occurs in the d-wave channel with a t-matrix given 
by 

£k,k',q(«w m ) = tq(iu m ) cos(2v5 k ) cos(2i/3 k /). (1) 
Here ioj m = 2mnT is the bosonic Matsubara frequency (un- 
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less noted otherwise, we take ft = 1, fcs = 1), and <pk = 
axctax\(ky/k x ). In the absence of an external magnetic field, 
the incoherent pair tunneling contribution I P (V) is given by 
4e times the imaginary part of the diagram shown in Fig. [j] (c) 

HI 



I p (V) = AeC 2 Sa 2 lmtq =0 (iu; m -> 2eV + iS), 



(2) 



where S is the junction area, a is the lattice spacing, and the 
coefficient C - which determines the magnitude of the pair 
current - is given by the following expression 



~N~ 2 



C = -^T^2 F p (iui n )G k (iu> n )G- k (-iu> n ) 

n , p , k 



x(|F klP | 2 ) imp cos(2^ k ) 



(3) 



Here we presume that the mechanism for electron transfer 
from A to B derives from impurity assisted hopping in the in- 
sulating layer separating A and B. We define m to be the num- 
ber of impurity scattering sites per unit area of the insulating 
layer, N is the number of sites of a layer, and (|Vp k | 2 )i mp is 
the impurity averaged single-electron transfer. The momenta 
p and k are two-dimensional vectors. In Eq. we have ne- 
glected the weak voltage dependence of C which is justified in 
the regimes we will be studying where the t-matrix dominates 
the voltage dependence of the pair current. 

To estimate the size of the pair current we have used the 
Bardeen-Cooper-Schrieffer (BCS) form for the Gor'kov func- 
tion to describe the superconducting A lead 



F p (iu> r , 



A A cos(2^ F 



e 2 p + A 2 A cos 2 (2^p) 



(4) 



Here A^ is the maximum of the associated d-wave supercon- 
ducting gap. The detailed nature of the Green's functions in 
B are of course very important in determining the particle- 
particle t-matrix. However, the coefficient C is obtained by 
summing over both the momentum and frequency variables 
of the propagators. Thus C is only marginally affected by the 
precise form of the single particle propagators. Whether one 
replaces the product of the B Green's functions by their non- 
interacting form 



G k (iw n )G_k(-iw„) 



1 



u 2 +el 



(5) 



or whether one uses an extreme limit [|17p of pseudogap theo- 
ries 



G k (iu; n )G- k (-iLL! n ) 



1 



"n + 4 + A l COS 2 (2</5 k ) 



(6) 



changes the estimate of C only by factors of order unity [118]. 
In these expressions e p and e k are the single particle energies 
in A and B, respectively. 

We will assume that the insulating layer gives rise to a dif- 
fuse [ 19 ,20] electron transfer with 



(I Vp k | 2 ) imD = \V \ 2 + |F 1 | 2 cos(2^ p )cos(2 ¥;k ) 



(7) 



More generally one could imagine expanding the impurity av- 
eraged single-electron transfer (| Vp k | 2 }i mp in two dimensional 
crystal harmonics. In Eq.(^J) we have kept only the uniform 
and d-wave pair transfer parts [|T]| . It is the second term in 
Eq.(Q) that will enter in our calculations. The required size 
of V\ will be discussed below, together with other junction 
requirements. Using Eqs. one finds that 

C = T?n i N A {n)N B {K)\V^A A TY J lA{u n )I B {u n ) , (8) 



where N A (0) and Ng(0) are the single-particle density of 
states per spin, per site for layer A and B, respectively, and 



cos 2 (2(^ p ) 
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A 2 A B cos 2 2(^ p 



■7T X /U)l+A AB 
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A,B 



(E{k n ) - K{k n )) 



(9) 



In the above expression k 2 = A 2 A b /(ijJ 2 + A 2 A B ), and K 
and E are the complete elliptic integrals of the first and second 
kind, respectively. Carrying out the Matsubara sum in Eq. (|^), 
we find that to within numerical factors of order unity 



C~ T n l 7V A (0)iV B (0)|y 1 | 



(10) 



Now at low temperatures, where A and B are both supercon- 
ducting, a similar calculation shows that the Josephson critical 
current is given by 



L = 2eC'A R S 



(ID 



with A b the low temperature maximum gap in B and the coef- 
ficient C is closely related to C given by Eq. ^ p2". Using 
this to normalize the strength of I p (V) we have 



I P (V) ^ Ej 



E, 



\mt{2eV) 



(12) 



Here Ej = ftl c /2e is the zero temperature Josephson cou- 
pling energy between A and B, E c = (S/a 2 )N B (0)A 2 B /2 is 
the condensation energy of the B cuprate layer, and Irm^u;) = 
NB{0)lmto(o-') is a dimensionless form of the t-matrix for 
q = 0. 

It can be seen from Eq. ( |l2| ) that the important quantity 
measured in a pair tunneling experiments is Imt(ui). The form 
of this function varies depending on the particular scenario 
adopted for describing the pseudogap. For a wide class of 
theories, Imi q (tj) can be expressed in terms of the pair sus- 
ceptibility Xq( w ) an d tne pairing coupling constant g 



Imtq(aj) 



.9 2 ImXq(w) 



[l+.gRe Xq M] 2 + [.9lmXqH 



(13) 
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A useful form of Eq. ( |13| ) for experimental comparison is dis- 
cussed in Refs. j7|,[§, although other alternatives may even- 
tually be proposed using different precursor scenarios 23 1 . 
The approach of Refs. provides a concrete diagrammatic 
prescription for computing \- For q = and sufficiently low 
frequencies, l+gRexo(w) « (a/j)(ui— oj ) and glmxo(u) ~ 
w/j, where (for T close to T C B ), w = ( 7 /a)(T/T c B - 1). 
The values of a and 7 depend on g. Under these conditions, 



Eq. (13) yields 



Imi(w) 



7W 



a 2 (u — lo ) 2 + LO 2 



(14) 



In the weak-coupling limit, where the dimensionless parame- 
ter a ~ 0, Eq. @ yields the well known result [|b[|T5|] 



w/7 



l) 2 + (c/7) 2 



(15) 



In this regime the pairing fluctuations are associated with crit- 
ical behavior and are essentially diffusive in nature. This case 
was addressed in earlier incoherent pair tunneling experiments 
Jl5| ] on conventional superconductors. By contrast, in the in- 
termediate coupling regime, which corresponds to a ~ 1, the 
value of u! a is strongly reduced, resulting in a pronounced res- 
onance in Imt(w) at this frequency. Thus, the pair fluctuations 
acquire a propagating nature [0|8]]. 

These two theoretical limits are illustrated in Fig. || which 
presents the self-consistently calculated t-matrix [^|] in weak 
and intermediate coupling, corresponding to B' and B respec- 
tively. Here T/T c s = 1.1. Notice the asymmetry InrF(w) ^ 
\mt{— uS) in the second case which provides a strong signature 
for pair resonance effects. This asymmetry is, in turn, related 
to an asymmetric density of states [Q], which may be associ- 
ated with that observed in STS experiments [[To|]. This figure 
also reflects the predicted voltage dependence of the pair tun- 
neling current. Within the superconducting pairing fluctua- 
tion scenarios of the type discussed in Refs. a prominent 
peak in ImJ(w) is expected to persist in underdoped cuprates 
to temperatures of order T ~ T*, considerably higher than 
T C B . Alternative scenarios [||^3|] can be used, presumably, to 
provide analogous signatures, within their respective theoret- 
ical framework. The importance of the incoherent pair tun- 
neling experiment lies in its ability to detect such features and 
therefore to confirm or falsify different classes of pseudogap 
scenarios. 

Let us now estimate the size of the pair current I p given 
by Eq.(|l2|). The condensation energy density can be in- 
ferred from heat capacity measurements For an under- 
doped Y Ba-zCuzO^+x of T c ~ 60K we obtain e con d ~ 
2 x 10 4 J/m 3 . Using typical values for junction surface area 



lQ-^m 2 [p 



and taking the layer thickness of order 

13 



s 

4 ~ 10 A, we find E c = e corid S£ c - 2 x 10 -13 J 
10 6 eV. Given the relatively large condensation energy, the 
strength of the Josephson coupling becomes crucial for the ef- 
fect to be observable. For a typical critical current of order 
I c ~ 10 mA, the corresponding Josephson coupling energy 



is Ej = hl c /2e - 3 x 10~ 18 J w 20 eV. Consequently, 
Ip ~ Ic(Ej/Ec) ~ 0.2 /j,A, which is of the same order as the 
pair currents detected in conventional superconductors [ |l5| ] . 

Thus it is important to fabricate junctions with c-axis 
Josephson current density in the range of 10 2 A/cm . Crit- 
ical current densities sustained by recently fabricated trilayer 
junctions [|l4j] are in this range. In order to detect this small in- 
coherent pair tunneling contribution to the total current I(V), 
one must be able to separate it from the larger quasiparticle 
current [j25|. Fortunately this can be done by turning on a 
magnetic field H in the plane of the junction, as shown in 
Fig. [I] (a). When H is such as to put several flux units in the 
junction, the incoherent pair tunneling will be suppressed. If 
the thickness of the A layer is larger than the in-plane penetra- 
tion depth, A„ b (T), and the B layer is thin compared to this, 
then one flux unit will be present when HLA Q b(T) « hc/(2e). 
For L = 100 /im and A^ b ~ 0.2 /.tm and d ~ 500A, one has 
H ~ 0.1 Gauss. In order to observe the zero field pair current 
I P {V), it will be necessary to magnetically shield the junc- 
tion. Then by subtracting the I-V data in the presence of an 
external magnetic field H « 1 Gauss from the I-V data in 
the absence of H one can obtain the incoherent pair tunneling 
contribution I P (V). Further complications might be caused 
by thermal voltage noise [ |l5| , [26| ] in the junction circuit due to 
the relatively elevated temperatures at which these measure- 
ments need to be carried out. One possibility would be to use 
the single-layer Bi2201 compound for both leads: this mate- 
rial has a phase diagram similar to that in Fig. [I] (b), but with 
a relatively low optimal T c . 

In conclusion, we have argued that the measurement of the 
pair tunneling current between an optimally doped and an un- 
derdoped cuprate can be used to probe the pairing fluctuations 
in the pseudogap state. This experiment has, in principle, the 
potential to reveal whether the pseudogap state is in fact due 
to pairing fluctuations. Indeed, strong pairing correlations in 
the pseudogap state will be manifest in a large pair current, 
as compared to the pair current of a junction where an over- 
doped lead of the same T c is used. No such strong doping 
dependence of the pair current is expected within pseudogap 
scenarios that do not invoke the onset of strong pairing corre- 
lations below T* . To illustrate this experiment, we have cho- 
sen the particular case of a c-axis junction geometry and iden- 
tified the region of the phase diagram where the experiment 
should be performed. We have also established the range of 
basic junction parameters suitable for observing the pair cur- 
rent, and suggested an experimental procedure for separating 
the small pair contribution from the large quasiparticle cur- 
rent. 
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FIG. 1. (a) Proposed experimental configuration for a junction 
involving two cuprate leads A and B, with transition temperatures 
as indicated by the phase diagram of (b). (c) Diagrammatic repre- 
sentation of the incoherent pair tunneling current contribution. Lines 
(±fc, ±fc') and ( ±p, ±p') correspond to single electron propagators 
of the normal pseudogapped (PG) B lead and anomalous Gor'kov 
propagators of the superconducting A lead, respectively. The dots 
(•) represent tunneling matrix elements Vpk, and the box stands for 
the particle-particle t-matrix of B. 




FIG. 2. Predicted voltage dependence of the pair tunneling cur- 
rent; following [^] the solid and dashed curves correspond to under 
(B) - and overdoped (B') leads at T/T c = 1.1; the slightly doping 
dependent Q, is of order of several hundred meV. Dashed curve is 
similar to that of conventional fluctuation picture Q13|. The asym- 
metry of the solid curve is an important signature which should be 
noted. 
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